Continuous Measurement Enhanced Self- Trapping of Degenerate Ultra-Cold Atoms in 

a Double- Well: Nonlinear Quantum Zeno Effect 
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In the present paper we investigate the influence of measurements on the quantum dynamics of 
degenerate Bose atoms gases in a symmetric double-well. We show that continuous measurements 
enhance asymmetry on the density distribution of the atoms and broaden the parameter regime 
for self-trapping. We term this phenomenon as nonlinear quantum Zeno effect in analog to the 
celebrated Zeno effect in a linear quantum system. Under discontinuous measurements, the self- 
trapping due to the atomic interaction in the degenerate bosons is shown to be destroyed completely. 
Underlying physics is revealed and possible experimental realization is discussed. 
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I. INTRODUCTION 

Double-well system is a paradigm to demonstrate 
quantum tunnelling phenomena and attracts much atten- 
tion from the societies of both theoreticians and experi- 
mentalists since the establishment of quantum mechanics 
For a symmetric double- well, the amplitude distribu- 
tions of all eigenstates are symmetric, so a particle is ex- 
pected to oscillate between the two wells and its chance to 
visit each well is even Q . However, there exist two effects 
that can break the above symmetry leading to an asym- 
metry on the particle's densitydistribution, i.e., quantum 
measurements [1, 0, S S S 11, H, fioft and collisions be- 
tween particles (many-body effect) [ill . [l2l Il3l . [lil [l5j . In 
the former case, the quantum measurement couples the 
system to outer environment and induces the quantum 
de-coherence of the system. As a result, the quantum 
tunnelling between two wells is suppressed completely 
so that the particle keeps staying in one well and has 
no chance to visit the other well. This is the celebrated 
quantum Zeno effect (QZE) @,G3. 

In the latter case, the 
two-body collisions among the degenerate Bose-Einstein 
condensates (BECs) atoms leads to a nonlinear excita- 
tion, manifesting a highly asymmetric density profile of 
the BECs even in a symmetric double-well. This some- 
how counterintuitive phe nomenon has been observed re- 
cently in labs [II, Hill. 

In the present paper, what we concern is how a mea- 
surement affects the dynamics of the many-body quan- 
tum system characterized by degeneracy and dilutencss. 
Under a mean-field approximation and without measure- 
ment, the dynamics of this system is described by a non- 
linear Schrodinger equation, known as Gross-Pitaevskii 
equation (GPE), where the nonlinearity arises from the 
interaction between the degenerated atoms and is respon- 
sible for the unusual self-trapping phenomenon. Consid- 
ering the measurement carried out by a position meter, 
the GPE is modified to be a stochastic nonlinear differ- 
ential equation. With solving the stochastic equation we 
achieve insight into how the measurement affects the self- 



trapping of BECs. Our main result is that continuous po- 
sition measurements enhance asymmetry on the distribu- 
tion of the atoms density profile and broaden the param- 
eter regime for self-trapping. We term this phenomenon 
as nonlinear quantum Zeno effect. However under dis- 
continuous measurements, the self-trapping phenomenon 
could be destroyed. Physics behind the above phenom- 
ena is revealed and possible experimental realization is 
discussed. 

This paper is organized as follows: We first briefly in- 
troduce the quantum measurement theory in Sec. II. In 
the third section, we describe our quasi 1-D double-well 
model of BECs. In Sec. IV, we discuss the influences of 
the continuous or discontinuous measurements on the dy- 
namics of the system. Finally, we summarize and discuss 
our work in section V. 



II. MEASUREMENT THEORY 

In the physics community, the measurement side of 
quantum theory is one of the fundamental issues 
Qj, [H. There exist various theories about this topic 
depending on concrete physical systems and measured 
physical quantities [HI, [2fj |. For the BECs system we 
consider, we make use of continuous (in time) posi- 
tion measurement theory [^, [1(3]. Within the frame- 
work of this theory, the position and the momentum of 
the system and the meter are denoted by {x(t),p(t}} 
and {X(t),P(t)}, respectively. Usually we can de- 
scribe the pseudo-classical meter as (Q\X (t) , P (t)) = 

(2^ 2 )~ 1/2 exp {- [Q ~ 4 X J t)]2 } exp {j-QP (t)} in the Q 

(position of the meter) representation, where a denotes 
the uncertainty of the meter pointer. Considering that 
the meter is pseudo-classical, one can treat this state as 
\X (f),0) 



\X(t),0) = (2ira 2 



-1/2 



exp 



[Q-x(t)f 

4ct 2 



(1) 
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This is one of the features of this theory that the higher 
up the chain from the system to the observer the cut is 
placed, the more accurate the model of the measurement 
will be. The interaction between the system and the 
meter can be described as 



H S m =jP[x- (x) s ] . 



(2) 



where (x) s is the average position of the system and 7 
represents the relaxation rate of the meter pointer to the 
system mean position (x) s . This interaction translates 
the position of the meter by an amount proportional to 
the average position of the system. The state of the 
system-meter at time t is taken to be 

|$(t)) = |*(*),o>®|*(t)>, 

where \$f (t)) is a system state vector. The combined 
state at time t + r is 

|$(t + r ))=expj-^P |X(t),0>®|* (* + 

(3) 



where 



!*,,(/ f -I) <:rp [ ) !*(,(/»> 



(4) 



III. DOUBLE- WELL MODEL 

We will focus on the case of BECs trapped in a quasi- 
one dimensional symmetrical double-well. Along the lat- 
eral directions the BECs is tightly confined so that the 
Hamiltonian governing the dynamics of BECs reduces to 
the following ID form, 



Ho 



\<P_ 

2dx 2 



+ V (x) + r)\ip 0)| 



(0) 



after re-scaling the energy unit by -^j^ and the length 
unit by L (effective size of the double- well potential); m 
is the mass of the trapped atom and r\ is the nonlinear 
parameter that is proportional to the 1-D reduced s-wave 
scattering length between the degenerate atoms and the 
total number of the atoms. As the size of the two well 
potential realized in the experiment [l3| is around L ~ 
13/im and the mass of the alkali atom is around m ~ 
T )'lO~ 25 kg, the order of the time in our paper should be 



From Eq.([3]) and Eq.([3]), we can see that the evolution 
of the combined system-meter is purely unitary. The 
meter has undergone the desired translation because of 
the interaction with the system. Meanwhile the system 
has undergone its usual free evolution combined with a 
measurement back-action from the meter pointer. 

After considering the read-out, the evolution of the 
combined equations in the first order of r can be written 
as a stochastic differential equation 

j t \*(t)) = {- l -H Q - T -[x-{x) s ] 2 

+Vf£(t)[s-<a>.])|*(t)) (5) 

where T — ^ 2 t j^a 2 denotes the interaction strength be- 
tween the meter and the system. The noise term £ (i) 
indicates Gaussian noise of standard deviation 1 and can 
be modelled by white noise (£ (t) £ (*')) = 6 {t - t'). 

The above stochastic Schrodinger equation is equiva- 
lent to a stochastic master equation for the selective evo- 
lution of the conditioned density operator Q . It is con- 
ditioned on the entire history of the meter readout X(t). 
If we were only interested in the non-selective evolution 
of the system, then we would have to discard all knowl- 
edge of the evolution of the system. This is achieved 
in the usual manner of averaging over all possible meter 
readouts at all times t. In our case this simply amounts 
to averaging over the stochastic term in Eq.(5), which 
gives zero. The process of photon scattering on a BEC 
that results in population difference measurement can be 
modelled by the non-selective evolution of the density op- 
erator [2lj . In our following discussions, for generality, 
we include both the de-coherence term and the noise term 
modelling the nondestructive measurements on BECs. 
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200 time duration in our dimensionless unit corresponds 
to 0.2s, that is within the lifetime of the BECs under the 
present experimental conditions. So, our calculations are 
extended up to time moment of 200 in following. Mean- 
while, the energy scale is around khz. The double-well 
potential is expressed by 



V(x) 
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To investigate the tunnelling dynamics of BECs in this 
double-well, we start from a superposition of ground 
and the first excited state of stationary CPE, i.e., 
Hoipo,i( x ) = A i o,i'0o,i(a ; )) where ^0,1 0*0 are ground state 
and the first excited wave function respectively and the 
corresponding chemical potential are /Xo,i- The coeffi- 
cients of the superposition are chosen so that the BECs 
are localized in one well initially, i.e., t = 0, as shown in 
Figffl 
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|*(x,0)) = -^=(V>oO)±V>iO)). 



(8) 



We should mention that, in addition to the above sym- 
metric states having linear counterpart, the CPE al- 
lows eigenstates that are totally asymmetric, so called 
symmetry-breaking states, in the strong interaction case. 
These asymmetric eigenstates are crucial in understand- 
ing the dynamics of BECs and are the source of unusual 
self-trapping phenomenon [22|, Our initial state is 

chosen as a superposition of ground and the first excited 
states and has asymmetric property, but it has nothing 
to do with the above asymmetric eigenstates. And fur- 
ther, because the superposition principle breaks down in 
the nonlinear case [23j |. we cannot predict the temporal 
evolution of a superposition state from the superposition 
principle. 
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FIG. 1: (Color online) The top two pictures show ground state 
and the first excited state, respectively. The bottom two show 
the wave function and the density profile of the initial sate, 
respectively. 



the operator-splitting method to calculate the evolution 
of the wave function and then calculate average position 
through the formula (x(t)) = J ip(x,t)*xil>(x,t)dx. In 
Figl^l we show the dynamic evolution of the coherent 
superposition state (jSj) for different nonlinear interaction 
parameters 77 by plotting the average position (x(t)} with 
respect to time. It shows quite different dynamic behav- 
iors, depending on the values of the nonlinear interaction. 
For the weak interaction case, i.e., i] < r\ c ~ 0.143, the 
BECs demonstrate a coherent Josephson oscillation be- 
tween two wells, while for r\ > r) c , the BECs are trapped 
in one well, a phenomenon so called as the self-trapping. 
In our case, the population difference between the two 
wells of our initial state is quite large, so the self-trapping 
condition is easier to be satisfied compared to usual cases 

ana. 

IV. THE INFLUENCE OF THE 
MEASUREMENT 




FIG. 2: (Color online) The Josephson Oscillation and Self- 
trapping effect at r = 



In this section, we will investigate the influence of 
the measurements on the dynamics of the system. Ac- 
cording to the measurement theory, the continuous mea- 
surements is modelled by treating a sequence of n mea- 
surement operations described in the Sec. I, each sepa- 
rated by a time interval At, over a total time duration 
T = nAt. The continuous limit is obtained by taking 
the limit At — > 9]. The probability that the sys- 
tem is found on the initial state at each independent 
measurement is Po(t = T) w 1 — 2(Aif) 2 -^-, where 
(AH) 2 = (ipo\H 2 \il)o) - (ipo\H\ifj ) 2 ■ As is shown in fol- 
lowing, one can recover the QZE in the linear dynamics 
of the system (77 = in © ) by setting the time step 
as At = 10~ 3 , which corresponds to the separated time 
r between two consecutive measurements. But if we in- 
crease the time interval to At = 10~ 2 , the continuous 
measurements condition is broken, we fall into the regime 
of the discontinuous measurements. 



In the absence of the nonlinear interaction between 
atoms, i.e., rj — 0, our Hamiltonian ([6]) reduces to one de- 
picting the motion of single particle in double-well poten- 
tial and the QZE is observed in the presence of continuous 
measurements [t| [To[ • When the nonlinear interaction is 
present and larger than a critical value r\ c 0, [HI, [r3 |. 
the somehow counterintuitive self-trapping phenomenon 
occurs even without performing measurement. In the fol- 
lowing part of this section, the above interesting behavior 
will be demonstrated. Under the combined nonlinear in- 
teraction and quantum measurements, the dynamics of 
BECs is dramatically influenced and detailed discussions 
will be presented in next section. 

Without performing measurements, i.e., T = 0, the 
dynamics of the system strongly relies on the interaction 
between the particles, the time-dependent GPE governs 
the evolution of the BECs. As in Ref. [12J, we apply 



A. Continuous measurements: The nonlinear Zeno 

effect 

With applying the continuous measurements, the dy- 
namics of the system is governed by Eqs.([5]). It is inter- 
esting to note that the equation ([5]) can be split into two 
parts. One is the normal free evolution part of the non- 
linear Schrodinger equation, the other part is the time- 
dependent random part. To solve the dynamic evolution 
of this problem, we first transform the wave function into 
the "interaction" representation and integrate the time- 
dependent part (the second part) with the Weiner incre- 
ment method. Then we apply the inverse operation to 
obtain the wave-function in the coordinate space. Please 
note that the transformation to the "interaction" rep- 
resentation is time-independent unitary transformation 
described as exp(— iHoAt/K), where At is the time step. 
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FIG. 3: (Color online) The average position of the BECs as 
the function of time with increase the value of T from the top 
to the bottom in the linear case. 



The Weiner increment dw helps us to numerically in- 
tegrate Eqs.© with a random term. According to the 
theory of Stochastic different equations, we introduce dw 
as follows, 



dw „ , . 



(9) 



Considering that our noise term is Gaussian white 
noise [24| , we can write Weiner increment as 



dw = \fr Rand (0, 1) 



(10) 



where Rand(0, 1) denotes the Gaussian random number 
with zero mean value and 1 standard variance. 

Besides the time-dependent random noise term, we 
also have to include the idea of the quantum trajec- 
tory theory [l^. In the quantum trajectory theory, 
each measurement operator can be regarded as a "state 
preparation" procedure. This procedure creates an en- 
semble of systems each possessing a certain state. Actu- 
ally, the measurement operator is realized by projecting 
the state into one of the position meter state ([1} . Physi- 
cally, we can understand this process as following: After 



finishing each measurement, we have destroyed the mea- 
sured quantum state and the state collapses into one of 
the states of meter pointer. So in the next step, we ex- 
actly start our evolution from this quantum state. There- 
fore, after each measurement quantum state is prepared 
through a normalization process. 

In Figj3l we plot the temporal evolution of average 
position of the particle for linear case rj = 0. It shows 
that, for weak measurement (i.e., small T), the parti- 
cle oscillates between the two wells, whereas the parti- 
cle will stay in the left well (its initial location) with 
increasing the interaction strength between the meter 
and the system. The above calculation clearly demon- 
strates the celebrated QZE. In the above calculation, the 
time step is set as At = 10~ 3 , which is small enough to 
guarantee that our measurements are continuous. With- 
out measurement, the period of the oscillation is around 
T os ~ 250. From our Eq.(5), the de-coherence term is 
estimated proportional to r/2 x (Axq) 2 , where Axo is 
the distance between two wells and is 2\^8. So it is ex- 
pected that when r~ (2x tt)/(16 x T os ) ~ 0.0015 the 
de-coherence effect caused by the measurement will be 
significant so that the linear QZE is observed. The above 
simple estimation excellently agrees with our numerical 
calculation as shown in Fig. 3, where we see at V > 0.001 
the particle turns to be localized in its initial well. 

However, we find that the above simple picture is no 
longer available when the nonlinear interaction between 
atoms emerges. In this situation, we find that the de- 
coherence effect caused by the measurement is signifi- 
cantly enhanced by the atomic interaction. 

We plot the average position of BECs in the presence 
of the nonlinear interaction between atoms, i.e., r\ = 0.1, 
in FigdJ In this case, the nonlinear interaction is smaller 
than the critical nonlinear interaction rj c , so that the 
BECs show a Josephson oscillation between the two wells 
in the absent of measurements T = (see Fig|2|). With 
increasing the interaction strength between system and 
meter, FigQ]shows that the BECs turn to stay in one well 
manifesting that the averaged position is smaller than 
zero. The critical value T for the occurrence of Zeno effect 
is modified from V = 0.001 in linear case to T = 0.0003 
in weak nonlinear case (77 = 0.1). In this case the period 
of the nonlinear Josephson oscillation is T os ~ 300. From 
the simple picture given in the above discussions we esti- 
mate that the de-coherence effect become significant only 
when T > 0.0012. However, at T = 0.0003 we already 
have observed the blockage of the Josephson oscillation. 
This analysis implies that the interaction among BECs 
atoms significantly enhance the de-coherence effect. As a 
result, even a moderately weak nonlinear interaction can 
play a dramatic role and significantly change the move- 
ment of the atoms. So we conclude that the nonlinear in- 
teraction broadens the regime of the localization of atoms 
and makes it easy to observe the quantum Zeno effect. 
Comparing with the linear quantum Zeno effect, we term 
this phenomenon as nonlinear quantum Zeno effect. On 
the other hand, the continuous measurements suppress 
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FIG. 4: (Color online) Same as Fig. ((3} in the case of r\ — 0.1 



FIG. 5: (Color online) Same as Fig[3]in the case of r\ c = 0.143 



the oscillation between the two wells and enhance the 
self-trapping by enlarging the regime of self-trapping. To 
show this point more clearly, we further plot the average 
position of the BECs with r\ c = 0.143 in Fig|5] 

When the nonlinear interaction is larger than the criti- 
cal value, we also calculate the average position of BECs 
with respect to time as shown in FigEl In the pres- 
ence of the strong nonlinear interaction, the BECs atoms 
are self-trapped even without performing the measure- 
ment. From Fig[6l we see that, the continuous measure- 
ments only cause small displace of the averages position 
of BECs, but can not drive the BECs to tunnel through 
the barrier. This means that the continuous measure- 
ments do not destroy this kind of nonlinear self-trapping. 

In FigJTl we plot the parameters diagram of the oc- 
currence of the self-trapping. We find that, the continu- 
ous measurements broaden the parameter regime for the 
emergence of self-trapping. For example, the continuous 
measurements of strength T = 0.0006 reduces the value 
of the critical nonlinear interaction rj c from 0.143 to 0.05. 
On the other hand, the nonlinear interaction between the 
degenerate atoms makes it easier to observe QZE by re- 
ducing the critical coupling strength T from 0.001 in lin- 



ear case (see Fig(3]) to 0.0003 in the case of r) = 0.1 (sec 
Fig|4j. Therefore, we conclude that, the combination of 
the continuous measurements and the nonlinear interac- 
tion greatly broaden the parameter regime for observing 
QZE and occurrence of the self-trapping. 



B. discontinuous measurement 

So far, we have investigated the effect of the contin- 
uous measurements on dynamics of the nonlinear quan- 
tum system. In this part, we want to know what happen 
if our measurement is not continuous, i.e., to study the 
influence of the discontinuous measurements on the self- 
trapping effect. 

To investigate the effect of discontinuous measure- 
ments, we will modify the time interval of the sequence 
of the measurement from 10~ 3 to 10~ 2 . Considering the 
barrier height, we can not increase the parameter T to too 
large. But if we keep r = and increase the nonlinear in- 
teraction to r\ = 0.5, the dynamics of BECs initially local- 
ized in one well demonstrates a self-trapping phenomenon 
the same that as shown in FigJ5] In FigJSJ we plot our re- 
sults in the same way as the previous figures from the top 



6 



g 
- < — * 

CO 

o 
D_ 

CD 
D) 
TO 
i_ 
CD 
> 
< 



-r=0.0002 
r=0.0003 



■r=0.0006 
r=0.0009 



-r=0.001 
r=0.005 



50 



100 

Time 



150 



200 



FIG. 6: (Color online) Same as Fig[3]in the case of r\ — 1.0 



to the bottom for various values of T. It shows that the 
oscillation amplitude of the averaged position of BECs in- 
creases with the increasing measurement strength T. In 
the case of T = 0.005, the particles can cross the barrier 
and tunnel into the other well. This means that the dis- 
continuous measurements break the self-trapping effect 
that is induced by the strong interaction between atoms. 
In the other words, the discontinuous measurements but 
still high frequency measurements have tendency to re- 
cover the macroscopic coherence properties of the system. 
The similar effect of the measurements process has been 
pointed out by investigating the dynamics of the quan- 
tum dynamics of the BECs with two wells including the 
nondestructive measurements 12511 . 



V. CONCLUSION 

In conclusions, we have investigated the influence of 
both continuous and discontinuous position measure- 
ments on the quantum dynamics of degenerate Bose 
atoms gases in a symmetric double- well. We find that, 
continuous position measurements enhance the asymmet- 
ric distribution of the atoms density profile and broaden 
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the parameter regime for self-trapping. Moreover, the 
de-coherence effect caused by the measurement is signif- 
icantly enhanced by the atomic interaction. The non- 
linear interaction between the degenerate atoms makes 
it easier to observe QZE by decreasing the critical cou- 
pling strength T, only over which observing QZE is possi- 
ble. Therefore, we conclude that, the combination of the 
continuous measurements and the nonlinear interaction 
greatly broaden the parameter regime for observing QZE 
and occurrence of the self-trapping. On the other hand, 
we find that, discontinuous measurements may break the 
self-trapping. It implies that the discontinuous position 
measurements may enhance the tunnelling of the BECs 
atoms. In the present experimental condition, BECs 
trapped in the symmetric double-well has been realized 
with using optical trap technique [l3| and nondestruc- 
tive measurements can be realized by shining a coher- 
ent light beam through BECs [25[. The measurement of 
the mean coordinate is, roughly speaking, equivalent to 
that of the inter-well population difference. Therefore, 
the photon scattering on BECs in a double-well poten- 
tial is sug gest ed to give the population difference mea- 
surement [2l| . Our predicted phenomena is hopefully 
observed with present experimental technique. 
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